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ABSTRACT 

Secondary  multiple-access  interference  processes  are  characterized  for  multi-hop 
packet  radio  networks,  in  which  users  are  assumed  to  be  Poisson-distributed  in  the  plane 
and  to  use  frequency-hopped  spread-spectrum  signaling  with  a  receiver-oriented  assign¬ 
ment  of  frequency-hopping  patterns.  The  throughput  per  node  and  the  average  forward 
progress  are  then  evaluated  for  frequency-hopped  multi-hop  networks  which  employ  (i) 
random  forward  routing  with  fixed  transmission  radius  (RFR)  and  (ii)  most  forward  pro¬ 
gress  routing  with  fixed  transmission  radius  (MFR).  The  optimal  average  number  of 
neighbors  and  transmission  radius  are  derived  for  these  cases  when  Reed-Solomon 
forward-error-control  coding  with  minimum  distance  decoding  or  binary  convolutional 
coding  with  Viterbi  decoding  is  employed. 


E.  Geraniotis  is  also  with  the  Systems  Research  Center  at  the  University  of  Maryland,  College  Park.  His  research 
was  supported  in  part  by  the  National  Science  Foundation  through  grant  ECS-85-16689  and  in  part  by  the  Naval 
Research  Laboratory. 

J.  W.  Gluck  is  an  NRL  fellow  under  ONR  Grant  N00014-85-G-0207. 


1 


I.  INTRODUCTION 

In  recent  years,  a  good  deal  of  research  has  been  done  in  the  areas  of  multi-user 
communication  networks  and  spread-spectrum  communications.  One  of  the  logical 
extensions  of  this  research  is  the  use  of  spread-spectrum  techniques  in  combination  with 
networking  techniques  in  order  to  provide  greater  multi-user  capabilities  and  higher 
resistance  to  interference,  whether  hostile  or  benign.  However,  there  have,  thus  far,  been 
few  attempts  to  pursue  this. 

The  general  feature  of  most  of  the  spread-spectrum  network  models  proposed  to 
date  has  been  a  lack  of  precision  in  the  characterization  of  the  effects  of  spread-spectrum 
techniques  on  network  performance.  The  general  practice  in  papers  like  [l]  and  [2]  has 
been  not  to  pay  sufficient  attention  to  the  accurate  modeling  of  the  effects  of  spread- 
spectrum  modulation,  therefore  the  results  obtained  can  be  very  optimistic  or  pessimis¬ 
tic,  depending  upon  what  assumptions  are  made.  In  [1],  a  model  is  presented  for  a 
frequency-hopped  (FH)  spread-spectrum  multiple  access  (SSMA)  digital  cellular  tele¬ 
phone  network,  and  an  expression  is  derived  for  throughput.  The  outstanding  features 
of  the  model  are  the  characterization  of  user  mobility  in  terms  of  a  two-dimensional 
Poisson  process  and  the  treatment  of  transmitter  power  attenuation  with  distance. 
Similar  techniques  are  also  used  in  [2],  in  which  the  authors  concentrate  on  trying  to 
model  spread-spectrum  multi-hop  networks  and  although  their  multi-hop  network  model 
has  merit,  their  treatment  of  the  spread-spectrum  modulation  is  inadequate  (they  only 
employ  the  processing  gain  without  taking  into  account  the  characteristics  of  the  specific 
spread-spectrum  modulation). 

In  contrast,  the  work  of  [3]  for  single-hop  frequency-hopped  spread-spectrum  net¬ 
works  examines  in  depth  the  effect  of  forward-error-control  coding  and  frequency-hopped 
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spread-spectrum  on  the  throughput  and  packet  error  probability  of  the  network  which 
employs  a  slotted  or  an  unslotted  ALOHA  protocol.  Furthermore,  the  work  of  [4]  inves¬ 
tigates  slotted  ALOHA  protocols  with  stable  throughput  for  frequency-hopped  Reed- 
Solomon  coded  single-hop  networks.  No  similar  attempts  to  analyze  accurately  spread- 
spectrum  multi-hop  networks  have  been  made  so  far. 

In  this  paper,  we  utilize  the  spread-spectrum  models  found  in  [3],  [4],  and  [8]  and 
incorporate  their  precise  detailed  characterizations  of  the  effects  of  frequency-hopping 
and  forward-error-control  coding  to  the  models  of  multi-hop  networks  treated  in  [2],  [5] 
and  [6].  The  effects  of  frequency-hopping  and  forward  error-control  coding  (Reed- 
Solomon  coding  with  minimum  distance  decoding  or  convolutional  coding  with  Viterbi 
decoding)  as  well  as  of  data  modulation  [M-ary  (or  binary)  orthogonal  FSK  with  non¬ 
coherent  demodulation]  are  taken  into  consideration,  and  the  interference  processes  are 
characterized  accurately  in  order  to  evaluate  the  throughput  per  node  and  the  average 
forward  progress  of  the  multi-hop  network. 

An  overview  of  the  contents  of  this  paper  is  as  follows.  In  Section  II,  results  are 
derived  for  the  throughput  per  node  and  the  average  forward  progress  for  frequency- 
hopped  multi-hop  networks  with  forward-error-control  employing  random  forward  rout¬ 
ing  with  fixed  transmission  radius  (RFR),  a  modification  to  the  scheme  treated  in  [2]. 
Section  III  contains  results  for  similar  multi-hop  networks  employing  most  forward  rout¬ 
ing  with  fixed  transmission  radius  (MFR),  adapted  from  [5],  Finally,  in  Section  IV,  we 


present  numerical  results  and  conclusions. 
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H.  FREQUENCY-HOPPED  MULTI-HOP  NETWORKS  EMPLOYING 
RANDOM  FORWARD  ROUTING  WITH  FIXED  TRANSMISSION  RADIUS  (RFR) 

In  this  section,  expressions  are  derived  for  throughput  and  expected  forward  pro¬ 
gress  (which  is  used  to  determine  optimum  transmission  ranges)  for  frequency-hopped 
multi-hop  networks. 

The  network  model  to  be  assumed  here  is  similar  to  those  found  in  [2],  [5]  and  [6], 

A  geographically-dispersed  network  is  assumed  with  a  uniform  traffic  matrix.  The 
number  of  users  in  a  given  area  is  a  two-dimensional  Poisson  process  with  density  X 
users  per  unit  area.  Each  node  is  assumed  to  be  able  to  transmit  to  users  located  within 
some  radius  R;  thus,  the  average  number  of  “neighbors”  (users  to  whom  a  given 
transmitter  can  transmit)  is  given  by  N  =  \nR2.  A  given  node  will  be  in  transmit 
mode  with  probability  p  and  will  not  be  in  transmit  mode  with  probability  1-p  (this 
represents  a  Bernoulli  process  with  parameter  p  ). 

The  channel  access  protocol  is  assumed  to  be  slotted  ALOHA,  and  the  system  is 
assumed  to  utilize  frequency-hopping  as  the  means  of  spreading  the  spectrum.  Receiver- 
oriented  assignment  of  frequency-hopping  patterns  is  assumed;  i.e.,  a  transmitter  with  a 
packet  to  send  to  some  receiver  uses  the  hopping  pattern  unique  to  that  receiver. 

Synchronization  at  the  packet  level  is  assumed  feasible  for  all  users.  Thus,  the 
uncertainties  in  the  timing  between  different  users  are  required  to  be  small  relative  to 
the  packet  duration;  however,  since  they  might  not  be  small  relative  to  the  dwell  time  of 
the  frequency-hopper,  we  consider  asynchronous  frequency-hopping  systems.  For  sys¬ 
tems  employing  Reed-Solomon  coding  one  codeword  per  packet  is  transmitted. 

There  are  many  possible  routing  strategies,  of  which  two  shall  be  analyzed  here: 
random  forward  routing  with  fixed  transmission  radius  (RFR)  analyzed  in  this  section 
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and  most  forward  routing  with  fixed  transmission  radius  (MFR)  analyzed  in  the  next 
section. 

Random  forward  routing  with  fixed  transmission  radius  (RFR)  is  based  on  the  com¬ 
pletely  random  routing  scheme  found  in  [2],  The  difference,  of  course,  is  that,  here,  a 
transmitter  always  chooses  a  receiver  in  the  forward  direction,  the  direction  in  which  the 
ultimate  destination  is  located. 

Figure  1  illustrates  this  strategy  and  the  types  of  interference  that  are  encountered 
therein.  In  the  figure,  X  is  transmitting  a  packet  that  is  destined  for  node  D.  X  ran¬ 
domly  chooses  a  node  Y  in  the  half-circle  of  radius  R  surrounding  X  that  is  in  the  direc¬ 
tion  of  node  D.  Forward  progress  is  denoted  by  the  length  of  segment 
XZ  —  z  =  rocos0o. 

In  the  figure,  X1  and  X2  interfere  with  X’s  transmission.  X1  is  a  primary 
interferer,  X1  has  a  packet  destined  for  D1,  and  he  randomly  chooses  Y  as  his  immedi¬ 
ate  destination  node,  the  same  node  chosen  by  X  (  note  that  X,  likewise,  is  a  primary 
interferer  with  Xx’s  transmission).  All  primary  interferers  use  the  same  spread-spectrum 
code  as  X  uses.  X2  is  a  secondary  interferer ;  X2  has  a  packet  destined  for  node  D 2  and 
chooses  Y2  as  an  intermediate  destination.  However,  since  the  distance  between  X2  and 
Y  is  less  than  R,  the  fixed  transmission  radius,  some  of  the  signal  power  used  by  X2  is 
received  by  Y .  The  spread-spectrum  codes  which  are  used  by  the  secondary  interferers 
are  different  from  the  code  used  by  X  for  transmitting  to  Y . 

In  analyzing  this  scheme,  what  is  desired  is  the  probability  that  X  transmits 
correctly  to  some  node  Y,  denoted  P(X—>Y).  In  order  to  evaluate  this  probability,  it 
is  necessary  to  evaluate  the  probabilities  of  occurence  of  a  number  of  contributing 


events. 
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Given  that  a  node  is  in  transmit  node  (with  probability  p ),  it  actually  transmits 
_N 

with  probability  1-e  2  ,  as  discussed  in  [3],  where  N  is  the  average  number  of  neighbor¬ 
ing  nodes  (within  a  circle  of  radius  R),  as  described  above.  Thus,  we  have 

PTx{X—*Y)  —  Pr  {X  in  transmit  mode}Pr  {X  can  find  a  receiver  |  X  in  transmit  mode} 
•Pr  (F  is  not  in  transmit  mode} 

This  accounts  for  the  “non-interference”  portion  of  P(X—*Y)\  we  now  concentrate  on 
the  interference-related  components. 

The  first  piece  of  information  necessary  in  order  to  evaluate  the  interference  com¬ 
ponents  is  the  number  of  potential  interferers  (surrounding  a  receiver).  The  probability 
that  there  are  k  interfering  users  (transmitting  packets)  within  a  radius  R  around  a 

e~N  Nk 

given  node  is  given  by  the  Poisson  distribution,  P  (K  =k )  = - - - . 

k ! 


(1-P). 


(1) 


Given  that  there  are  k  potential  interferers,  we  must  now  account  for  the  fact  that 
each  of  these  can  either  not  interfere  (i.e.,  not  transmit),  be  a  primary  interferes  or  be  a 
secondary  interferer.  Let  us  say  that  j  of  these  k  do  interfere;  then  k-j  do  not,  and 


this  probability  is  given  by 


n' 

1  -p 

l-e  2 

..  > 

1  *-/ 


Now,  say  that  i  out  of  the  j  are  pri¬ 


mary  interferers  and  j —  i  are  secondary  interferers.  The  probability  of  this  being  true 

V-i 


can  be  shown  to  be  X- 
N 


P  l-e 


n 

2  _P_ 
N 


The  next  step  involves  the  evaluation  of  the  probability  that  X  (out  of  the  *  +1 
users  transmitting  to  Y)  is  captured  by  the  receiver  Y  and  the  probability  that  Y 
correctly  receives  the  packet  of  X  in  the  presence  of  the  i  primary  and  j-i  secondary 
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interferes.  Regarding  capture,  we  assume  that  perfect  capture  occurs;  that  is,  only  one 
of  the  contending  users  will  be  captured  by  the  receiver,  and  this  event  has  probability 
1.  Since  i+l  users  (including  X)  are  contending  and  each  one  is  equally  likely  to  be 

captured  by  the  receiver  Y ,  the  probability  that  Y  captures  X  is  — .  Regarding  the 

t  +1 

probability  of  correct  reception  in  the  presence  of  j-i  secondary  interferers  it  will  be 
denoted  by  Pc  (j -i),  which,  for  various  coded  systems,  is  given  in  Appendix  A. 


Then,  summing  over  all  possibilities, 


00  e~N  Nk  k 
Pi  I  tx(X—*Y)  =  E  ~  7T~~  S 


“n  k! 


k=  0 


‘  £ 

,  =o  I  *1 


y=o  uj 


1-p  1-e 


N_ 

2 


1 


N 


N  i+i  ’  N 


}-* 


PAj-i)- 


(2) 


We  now  take  our  two  components,  Ptx(X—*Y)  and  Pj  |  TX(X  —+Y),  and  get  the 
overall  expression  for  P(X—>Y).  This  is  simply 


P(X-^Y)  -  Ptx(X-+Y)P >  |  tx(X-+Y).  (3) 

Unfortunately,  the  — - —  in  Pj  \  jX{X—*Y)  presents  an  obstacle  in  any  attempt  to 
i  +1  1 

write  a  closed-form  expression  for  P (X  —*Y).  However,  -t— —  can  be  closely  approxi- 

t  +1 

j  * 

mated  by  a  sum  of  exponentials,  J]  c„e  ",  in  which  the  Cj,’s  and  7„’s  can  be  determined 

V 

using  Prony’s  method  (see,  for  example,  [7,  pp.  378-382]).  We  have  used  a  three  term 
approximation  in  our  analysis  below  and  in  computing  the  numerical  results  shown  in 
Section  IV. 

If  we  further  assume  that  Pc  (j-i)  can  be  put  in  the  form  E  bm,i  am,l  (this  is  true, 

m,l 

for  example,  for  coded  systems,  as  discussed  in  Appendix  A),  then  P(X—*Y)  can  be 
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expressed  in  a  closed  form  as 


N 


P(X^Y)  =  p\  l-e 


(1  ~P) 


'  Yj  c  v  Yj  X; 

i^=l  m  I 


,1  expj- 


pN 


N 


N 


O  ~~  Q  1 

-  -N~a"J  ^  -N 


) 


(4) 


where  8V=  e  . 


In  this  case,  the  local  normalized  throughput  (in  packets  per  hop,  per  node,  fre- 


r  log 2M  i 

quency  slot  and  dimension)  is  just - — — P(X—*Y),  where  we  multiplied  by  — , 

q  M  q 


log  -M 
M 


and  the  forward-error-control  (FEC)  code  rate  r  that  is  used  to  account  for  the 


bandwidth  spread.  The  expected  forward  progress  Z  can  be  shown  to  be  given  by 


Z  —  —P  {X  -+Y).  (5) 

7T 

In  Appendix  A  expressions  for  P(X—*Y)  which  permit  the  evaluation  of  the 
troughput  per  node  and  the  average  forward  progress  are  derived  for  Reed-Solomon  cod¬ 
ing  with  error-correction  decoding  and  erasure/error-correction  minimum  distance  decod¬ 
ing  s  well  as  for  binary  convolutional  codes  with  side  or  without  side  information  and 
Viterbi  decoding. 


As  an  example  of  a  case  in  which  Pc  (j—i)  is  of  the  form  bm  /  J  which  results 

m  ,1 


in  a  closed  form  expression  for  for  P  (X  —*Y),  we  consider  the  case  of  RFR  with 
RS(n,k)  coding  used  along  with  error-correction  decoding.  From  Appendix  A  eq.  (A. 2) 
(see  also  [3]),  we  notice  that  we  can  put  Pc(j-i )  in  the  desired  form,  where 


*"m  ,1 


(1  -Ph)n+l  m  and  b 


/  (-!)'( 1-P0) 


Thus,  (2)  becomes 


8 


N 


P(X^Y)  =  p\l-e  2  (1  -p) 


3  t  m  n  m 

E  E  E  #  (-1)'  d-Po) 

J/=1  m  =0  /  =0  m  ‘ 


mc  (n  +l-m ) 


J{ 


exp\  -pN 


1-e 


.E  c 
2 

~  ~N 


(1  ~Pk) 


n  +1  -m 


N 


(6) 
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HI.  FREQUENCY-HOPPED  MULTI-HOP  NETWORKS  EMPLOYING 
MOST  FORWARD  ROUTING  WITH  FIXED  TRANSMISSION  RADIUS  (MFR) 

This  strategy  differs  from  RFR  in  that  the  node  Y  that  results  in  the  most  forward 
progress  is  always  chosen  by  X .  This  scheme  is  discussed  in  [5],  The  fact  that  causes 
the  analysis  of  MFR  to  differ  from  that  of  RFR  is  that,  as  shown  in  Figure  2  since  Y  is 
the  “most  forward”  node,  there  can  be  no  nodes  beyond  Y  in  the  direction  of  the  desti¬ 
nation;  thus,  a  region,  denoted  A  ( r0,60 )  is  excluded. 

In  [5]  expressions  are  derived  for  fT  ^(v)  and  fz(  ),  the  distributions  of  the  position 
of  the  immediate  destination  node  {Y ,  in  the  case  of  X  or  a  primary  interferer;  some 
Y 2,  in  the  case  of  a  secondary  interferer).  Expressions  for  the  area  of  A  ( r0,90 )  are  deter¬ 
mined,  as  well.  These  are  all  needed  in  the  evaluation  of  P{X—*Y). 

As  a  preliminary  to  the  analysis,  define  A'  ( r0,60 )  —  nR  2-A  (ro,0o)  to  be  the  area 
of  the  non-excluded  region.  Thus,  N—\nR2  is  the  expected  number  of  neighbors  in  an 
entire  circle  of  radius  R,  while  N'  ( r0,d0)—\A '  ( r0,60 )  would  be  the  expected  number  of 
neighbors  in  the  non-excluded  region. 

The  quantities  to  be  determined  in  this  analysis  are  basically  identical  to  those  in 
the  analysis  of  RFR,  but  there  are  some  significant  differences.  One  such  difference  is 
that  the  Poisson  distribution  uses  N'  ( ro,0o ),  rather  than  N ,  in  the  distribution  of  the 
number  of  neighbors  around  Y. 

In  evaluating  the  interference  probabilities  in  this  case,  we  will  separate  the  proba¬ 
bility  of  there  being  j  inteferers  out  of  k  potential  interferers,  given  by 

from  the  probabilities  of  there  being  j-i  secondary  interferers  and  t  primary  interferers. 
These  two  quantities  are  now  conditioned  on  there  being  j  interferers. 


P( l~e  2) 
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If  an  interferer,  say  X1,  is  a  primary  interferer,  then  Y  must  be  X1(s  “most  for¬ 
ward”,  node.  The  overall  probability  of  a  node  being  a  primary  interferer,  given  that  it  is 
an  interferer  (i.e.,  that  it  transmits  and  is  a  neighbor),  is  given  by 

P^R  )  =  Pr  { X 1  is  primary  |  X1  interferes} 

R 

—  J  P  {Y  is  the  most  forward  neighbor  of  X1  |  Y  is  at  radial  distance  rx  from  X1} 
P  {X1, s  radial  distance  from  Y  is  between  rx  and  r  ^dr  t}dr  1  (7) 


Thus,  we  can  write 


K  2r  j 
~R^ 


J0  ‘  J„  fr/r'  ,e'  )dr  do' 


dr  j 


(8) 


2r  i 

where  the  term  in  brackets  and  respectively,  denote  the  first  and  second  probabili¬ 

ty  terms  in  the  right  member  of  (7). 

The  probability  that  some  node,  say  X2,  is  a  secondary  interferer,  given  that  it  is 
an  interferer,  is  given  by  the  complement  of  (8);  i.e., 


Pr  (X2  is  secondary  |  X2  interferes}  =  P 2{R  )  =  1  -  P  X{R  )  (9) 

This  quantity  represents  the  sum  of  the  probabilities  that  Y  is  not  in  X2’s  forward 
direction  and  that  Y  is  in  X2’s  forward  direction  but  is  not  “most  forward”. 


As  before,  we  now  sum  over  all  possibilities,  and,  unlike  the  RFR  case,  we  must 
now  also  take  the  expectation  over  r0  and  60,  the  position  of  the  “most  forward”  node. 
This  gives  the  following  result: 

JL  »  \  ■„  e-K' ['V' <r  „,(?„)]* 

P(X-r)  =  p( l-e  !)(1  -P)/0  J,E  - lT, - — 
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P(  1-e  )  1-P  1-e 


[1  -P1(R  W^PcU-i) 


■fr/r0,d0)de0dr0. 


In  deriving  (10),  we  assumed  that  capture  occurs  with  probability  1,  that  is,  one  of  the 
t +1  contending  packets  is  captured  by  the  receiver.  Furthermore,  if  we  use  the  three 

term  exponential  approximation  to  — - —  and  assume  that  Pc(j-i)  is  of  the  form 

i  +1 

X)  ^ m,iam~l<  35  we  did  in  the  RFR  analysis,  the  integrand  can  be  simplified  to  some 

m.l 

extent,  by  evaluating  the  summations.  At  the  end  of  this  section,  we  give  an  example  in 
which  this  evaluation  can  be  carried  out  to  the  end. 

r  log2Af 

Again  the  normalized  local  throughput  is - — —  P(X-*Y )  and  the  expected  for- 

q  M 

ward  progress  can  be  evaluated  by  multiplying  the  integrand  in  (10)  by  rocos0o  and  then 
integrating  as  in  (10). 

As  an  example,  when  Reed-Solomon  coding  with  error-correction  minimum  distance 
decoding  is  employed  we  can  obtain  for  the  MFR  case  (see  Appendix  B): 


P{X-+Y)  =  p\  1-e  J  (1-p) 


■hc»h  e 

v—l  m—  0  1=0  V  /v  / 


+l-m ) 
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R  \ 

f  f  exp 

JQ  J  It 
~  2 


~pN  ( r0,60)\  1-e 


i -SfAR  Hi -P*) 


n  +/ -m 


[i-Pi(R )]] 


fr,e(r  oJo)ddOdrO- 

Expressions  similar  to  (11)  can  be  derived  for  the  case  of  erasure/error-correction 
decoding  of  Reed-Solomon  codes,  as  well  as  for  the  Viterbi  decoding  of  binary  convolu¬ 
tional  codes.  These  expressions  have  been  derived  and  cited  in  Appendices  B  for  the 
MFR  case. 
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IV.  NUMERICAL  RESULTS 


In  this  section,  results  are  presented  for  the  normalized  throughput  per  node  and 
the  expected  forward  progress  for  frequency-hopped  multi-hop  networks  using  the  RFR 
and  MFR  schemes,  discussed  in  Sections  II  and  III,  respectively. 


The  usual  performance  measures  in  multi-hop  networks  (e.g.,  see  [5]- [6])  are:  the 

r  log2M 


normalized  local  throughput  defined  as  5  = 


-P(X— >Y)  (in  packets  per  hop,  per 


q  M 

node,  frequency  slot  and  dimension),  the  expected  forward  progress  Z  (in  distance  units 
per  hop  where  Z  has  the  same  units  as  the  transmission  radius  R ),  and  the  normalized 
total  throughput  7  (in  packets  per  node,  frequency  slot  and  dimension),  which  is  defined 


as  7  =  Z  45n*/  ——■S.  Notice  that  the  total  throughput  can  be  easily  obtained  from 


the  local  throughput  and  the  expected  forward  progress;  we  do  not  present  any  results 
on  the  total  throughput  in  this  paper.  Also  notice  that  the  expected  forward  progress 
and  the  total  throughput  depend  on  both  X  and  N ;  in  contrast,  the  local  throughput 
depends  on  X  only  through  N  =  \nR2,  the  average  number  of  neighbors. 

We  start  the  presentation  of  results  with  the  RFR  case.  Figure  3  shows  the  depen- 
dance  of  Z  on  N  as  parametrized  by  X.  The  expected  forward  progress  first  increases 
and  then  decreases  as  N  increases  for  fixed  X.  As  X  increases  and  N  is  held  fixed,  the 


forward  progress  decreases  because  the  radius  of  transmission  R 


decreases. 


Figures  4(a)  and  4(b)  show  that  both  S  and  Z  are  affected  by  the  value  of  p  ;  the 
values  increase,  then  decrease  with  p  .  The  optimum  value  of  p  when  N  is  an  unres¬ 
tricted  parameter  has  been  computed  to  be  approximately  0.466  for  S  in  Figure  4(a)  and 
approximately  0.355  for  Z  in  Figure  4(b).  The  value  of  p  =  0.1  appears  to  provide 
almost  invariant  5  or  Z  independent  of  the  value  of  neighbors  N .  One  thing  should  be 
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noted  regarding  p  and  X:  these  parameters  reflect  users’  behavior  and  can  not,  in  gen¬ 
eral,  be  set  to  desired  values. 

Figures  5(a)  and  5(b)  demonstrate  the  effects  of  changing  the  rate  of  the  code  and 
the  decoding  algorithm.  As  shown,  throughput  decreases  significantly  when  the  code 
rate  decreases  while  it  remains  invariant  for  a  larger  range  of  N  when  erasure/error- 
correction  decoding  rather  than  error-correction  decoding  is  employed.  In  contrast  the 
expected  forward  progress  increases  drastically  when  the  code  rate  decreases  and  even 
more  so  when  erasure/error  decoding  instead  of  error-only  decoding  is  used.  The  decrease 
to  throughput  due  to  the  decrease  in  the  code  rate  is  due  mainly  to  the  normalization  of 
the  throughput;  that  is,  the  increase  in  the  probability  of  correct  reception  due  to  the 
lower  rate  does  not  balance  the  penalty  of  multiplying  the  throughput  by  the  code  rate. 
On  the  other  hand  the  expected  forward  progress  does  not  involve  any  normalization  so 
that  it  takes  advantage  of  the  full  benefit  of  lower  rate  coding.  In  comparing  two 
schemes  with  similar  parameters  but  error-correction  capability  (e.g.,  different  code  rates 
or  decoding  schemes)  maintaining  a  large  value  of  throughput  over  a  broader  range  is  a 
manifestation  of  the  fact  that  the  one  scheme  rejects  the  multiple-access  interference 
better  than  the  other. 

Figure  6  shows  the  behavior  of  the  throughput  S  as  a  fuction  of  N  when 
parametrized  by  X  for  frequency-hopped  multi-hop  networks  employing  the  binary  con¬ 
volutional  code  of  constrained  length  7  and  rate  1/2  with  no  side  information  available 
to  the  decoder.  Similar  observations  as  for  Figure  3  are  in  order  here.  Notice  that  the 
values  of  Z  in  Figure  6  are  smaller  than  the  corresponding  values  of  Figure  3  and  for 
fixed  X  they  remain  invariant  with  respect  to  N  in  a  smaller  range. 
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Figures  7(a)  and  7(b)  show  S  and  Z,  respectively,  as  functions  of  N  when 
parametrized  by  p  for  the  same  binary  convolutional  code  as  above.  In  particular,  Fig¬ 
ure  7(a)  shows  the  same  trends  as  Figure  4(a)  does;  the  difference  lies  in  that  the  convo¬ 
lutional  code  provides  larger  peak  throughput  (for  fixed  p )  but  it  maintains  it  for  a 
smaller  range  of  N  than  the  Reed-Solomon  code.  In  contrast,  Figure  7(b)  shows  that 
the  convolutional  code  both  provides  a  smaller  peak  expected  forward  progress  and  it 
maintains  for  a  smaller  range  of  N  than  the  Reed-Solomon  code. 

Figures  8(a)  and  8(b)  show  S  and  Z ,  respectively,  as  functions  of  N  for  binary 
convolutional  codes  of  rates  1/2  and  1/3  and  constraint  lengths  7  and  9.  We  observe 
that  decreasing  the  rate  of  the  code  lowers  the  throughput  and  increases  drastically  the 
average  expected  progress.  Increasing  the  constraint  length  improves  the  performance  of 
the  network  in  terms  of  both  the  throughput  and  the  average  expected  forward  progress. 

Figures  9(a)  and  9(b)  provide  a  comparison  of  S  and  Z ,  respectively,  as  functions 
of  N  for  the  rate  1/2  binary  convolutional  codes  of  constraint  lengths  7  and  9  and  the 
Reed-Solomon  code  (32,16);  binary  FSK  modulation  with  noncoherent  demodulation  and 
hard  decisions  are  employed.  Figure  9(a)  shows  that  although  the  peak  throughput  of 
the  convolutional-coded  network  is  slightly  larger  than  that  of  the  Reed-Solomon-coded 
network,  the  Reed-Solomon-coded  network  maintains  a  larger  value  of  the  throughput 
for  a  broader  range  of  values  of  N .  Then,  as  Figure  9(b)  shows  the  RS-coded  network  is 
far  superior  to  the  corresponding  convolutional-coded  networks  in  terms  of  both  the 
peak  value  and  the  range  of  N  for  which  the  peak  value  is  approximately  maintained  of 
the  expected  forward  progress  Z . 

We  now  move  on  to  discuss  the  results  for  the  MFR  case.  Figures  10(a)  and  10(b) 
are  analogous  to  Figures  4(a)  and  4(b);  they  show  the  throughput  and  expected  forward 
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progress,  respectively,  parameterized  by  p  .  The  same  trends  that  appear  in  the  RFR 
case  also  appear  here  in  the  MFR  case. 

Figures  11(a)  and  11(b)  show  the  throughput  and  the  expected  forward  progress  for 
two  Reed-Solomon-coded  networks  with  code  rates  1/2  and  1/4.  The  throughput  of  the 
RS(32,16)-coded  network  is  higher  than  that  of  the  RS(32,8)-coded  system  for  most 
values  of  N  this  is  reversed  for  N  larger  than  80.  On  the  other  hand  the  expected  for¬ 
ward  progress  of  the  lower  code-rate  scheme  is  higher  than  that  of  the  higher  code-rate 
scheme  for  all  values  of  N . 

Finally,  Figures  12(a)  and  12(b)  compare  RFR  and  MFR  performance  in  terms  of 
throughput  and  average  expected  progress,  respectively,  for  a  given  set  of  parameters. 
In  Figure  12(a),  it  is  shown  that  the  RFR  scheme  with  RS(32,16)  and  error-only  decoding 
outperforms  in  terms  of  peak  throughput  value  the  corresponding  MFR  schemes  with 
the  same  RS  coding  and  error-only  or  erasure/error  decoding;  however  the  MFR  schemes 
maintain  a  large  value  of  the  throughput  for  broader  range  of  values  of  N  and  between 
the  two  the  erasure/error  decoding  scheme  is  superior.  Then  in  Figure  12(b)  it  is  shown 
that  the  MFR  scheme  outperforms  the  corresponding  RFR  scheme  and  that  for  MFR 
schemes  the  value  of  Z  can  improve  considerably  by  using  erasure/error  decoding. 
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V.  CONCLUSIONS 

We  have  presented  here  an  analysis  of  a  model  for  randomly-distributed  frequency- 
hopped  multi-hop  networks.  The  aspect  of  this  analysis  which  differs  most  from  previ¬ 
ous  analyses  of  multi-hop  spread-spectrum  networks  is  the  characterization  of  other-user 
interference  and  its  mitigation  by  the  use  of  spread-spectrum  techniques  (frequency¬ 
hopping  here).  We  have  also  been  able  to  incorporate  the  effects  of  modulation,  coding, 
and  noise  into  the  analysis. 

The  numerical  results  presented  show  network  performance  (normalized  local 
throughput  and  expected  forward  progress)  as  a  function  of  the  average  number  of 
neighbors  around  a  node.  Variations  in  performance  are  demonstrated  for  variatins  in 
probability  of  a  user  being  in  transmit  mode,  density  of  users  in  the  plane,  and  coding 
employed.  These  results  also  compare  the  two  types  of  routing,  RFR  and  MFR,  con¬ 
sidered  in  the  analysis. 

An  obvious  extension  of  this  work  would  be  to  consider  deterministic,  rather  than 
randomly-distributed,  networks.  This  would  not  involve  a  major  modification  to  the 
analysis;  it  would  involve  removing  the  Poisson-related  elements  in  the  analysis  and  is, 
thus,  rather  trivial.  However,  it  would  be  somewhat  more  difficult  were  one  to  consider 
non-identical  nodes,  i.e.,  nodes  with  different  transmitter  powers  (or,  similarly,  to 
account  for  signal  attenuation  with  distance). 

The  importance  of  this  work  should  not  be  underemphasized;  however,  there 
remains  much  to  be  done  in  the  area  of  spread-spectrum  multi-hop  networks.  The 
analysis  here  does  not  deal  with  issues  like  spreading  code  protocols  and  access  protocol 
stability.  With  so  many  unsolved  problems,  this  remains  a  fertile  area  of  research. 
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Figure  4(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  p  (RFR,  A=0.1,  RS(32,16) 
coding  with  error  -only  correction  decoding,  q=100,  N^=10,  SNR=20  dB,  AWGN,  32-ary 
FSK  with  non-coherent  demodulation,  asynchronous) . 
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Figure  4(b).  Expected  Forward  Progress  versus  Average  Number  of  Neighbors  for  Various  p  (RFR, 

A=0.1,  RS(32,16)  coding  with  error  -only  correction  decoding,  q=100,  N^=10,  SNR=20  dB 
AWGN,  32-ary  FSK  with  non-coherent  demodulation,  asynchronous). 


2.0x10' 


J — L-l 


Figure  5(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  Reed-Solomon  Code  Rates 
and  Decoding  Schemes  (RFR,  p=0.3,  A=0.1,  q=100,  N  =10,  SNR=20  dB,  AWGN ,  32-ary  FSK 
with  non— coherent  demodulation,  asynchronous). 
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Figure  5(b).  Expected  Forward  Progress  versus  Average  Number  of  Neighbors  for  Various  Reed-Solomon 
Code  Rates  and  Decoding  Schemes  (RFR,  p=0.3,  A=0.1,  q=100,  N  =10,  SNR=20  dB,  AWGN, 
32-ary  FSK  with  non-coherent  demodulation,  asynchronous).  b 


Figure  6.  Expected  Forward  Progress  versus  Average  Number  of  Neighbors  for  Various  A  (RFR, 
p=0.1,  CC(7,l/2)  coding  with  no  side  information  available  to  the  decoder,  q=100 
N  =10,  SNR=20  dB,  AWGN,  binary  FSK  with  non-coherent  demodulation,  asychronous) . 


Figure  7(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  p  (.RFK,  A=u.l, 

coding  with  no  side  information  available  to  the  decoder,  q=100,  N^=10,  SNR=20  dB 
AWGN ,  binary  FSK  with  non-coherent  demodulation,  asynchronous). 
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Figure  8(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  Binary  Convolutional 

Codes  (p=0.3,  X=0 . 1 ,  no  side  information  available  to  the  decoder,  q=100,  N^=10, 
SNR=20  dB,  AWGN ,  binary  FSK  with  non-coherent  demodulation,  asynchronous). 
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Figure  8(b).  Expected  Forward  Progress  versus  Number  of  Neighbors  for  Various  Binary  Convolutional 

Codes  (p=0.3,  A=0.1,  no  side  information  available  to  the  decoder,  q=100,  N^=10,  SNR=20  dB, 
AWGN,  binary  FSK  with  non-coherent  demodulation,  asynchronous) . 
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Figure  9(b).  Expected  Forward  Progress  versus  Average  Number  of  Neighbors:  Comparison  of  Binary 
Coding  Schemes  (RFR,  p=0.3,  A=0.1,  no  side  information  available  to  the  decoder 
(error  —only  correction  £n  the  Reed-Solomon  case),  q=.lOO,  N^=10,  SNR— 20  dB,  AWGN, 
binary  FSK  with  non-coherent  demodulation,  asynchronous). 


Figure  10(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  p  (MFR ,  1=0.1,  RS(32,16) 
coding  with  error  -only  decoding,  q=100,  N  =10,  SNR=20  dB,  AWGN,  32-ary  FSK  with 
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Figure  10(b).  Expected  Forward  Progress  versus  Average  Number  of  Neighbors  for  Various  p  (MFR, 
A=0.1,  RS(32,16)  coding  with  error  -only  decoding,  q=100,  N^=10,  SNR=20  dB,  AWGN, 
32-ary  FSK  with  non-coherent  demodulation,  asynchronous). 


Figure  11(a).  Throughput  versus  Average  Number  of  Neighbors  for  Various  Reed-Solomon  Code  Rates 
(MFR,  p=0.3,  X=0.1,  error  -only  decoding,  q=100,  N^=10,  SNR=20  dB,  AWGN,  32-ary 
FSK  with  non-coherent  demodulation,  asynchronous). 


Figure  11(b).  Expected  Forward  Progress  versus  Average  Number  of  Neighbors  for  Various  Reed-Solomon 
Code  Rates  (MFR,  p=0.3,  1=0.1,  error  -only  decoding,  q=100,  N^=10,  SNR=20  dB,  AWGN, 
32-ary  FSK  with  non-coherent  demodulation,  asynchronous) . 
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Figure  12(a).  Throughput  versus  Average  Number  of  Neighbors  for  Comparable  RFR  and  MFR  Systems 
(p=0.5,  X=0 . 1 ,  q=100 ,  N  =10,  SNR=20  dB,  AWGN ,  32-ary  FSK  with  non-coherent 
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Appendix  A:  Derivation  of  Expressions  for  Coded  RFR  Systems 

In  Section  II  a  result  is  stated  for  a  Reed-Solomon  coded  RFR  system.  In  this 
appendix,  it  will  be  shown  exactly  how  this,  as  well  as  other  results  for  coded  RFR  sys¬ 
tems,  can  be  determined. 

In  Section  II  a  general  result  for  RFR  systems  is  derived.  We  have  that 


P(X-*Y)  =  Ptx{X-+Y)Pi/tx{X-+Y) 


=  pi  l-«  (1-p)  E 
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(A.l) 


In  that  same  section,  it  is  stated  that  - -  can  be  closely  approximated  by  a  sum  of 


i  +1 


exponentials  of  the  form  J6„,  where  6„  —  t  ‘'in  the  notation  used  in  that  section.  A 

V 

three  term  approximation  is  used  here. 


In  order  to  obtain  a  closed-form  expression,  we  must  have  Pc{j-i)  such  that  the 

dependency  on  j-i  is  exponential;  for  example,  if  Pc  (j-i )  =  E^m ,/  am~l  >  this  require- 

m  ,1 

ment  would  be  satisfied.  Having  Pc{j-i)  in  such  a  form,  results  in  the  final  sum  in 
(A.l)  having  the  form  of  a  binomial  expansion.  Several  coding  schemes,  including  those 
used  here  (Reed-Solomon  and  binary  convolutional),  result  in  a  Pc(j-i),  of  the  required 
form. 

We  shall  first  consider  Reed-Solomon  coding  RS (n  ,k )  with  error-correction 


decoding.  In  this  case,  we  have 


2 


Pc(j-i)= 


(A.2) 


where  t  = 


is  the  error-correction  capability  of  the  RS(n  ,k)  code  and  pe{j-i)  is 


the  probability  of  a  symbol  error.  Furthermore,  ps(j-i)  can  be  bounded/  approximated 
by 


PBU-i)  <  1-(1  -P0)(l-Phy  ‘  (A. 3) 

for  codes  over  GF(Ml)  [(1-P0)  becomes  (1  -P0)m‘  for  codes  over  GF  where  P0  is 

the  probability  of  symbol  error  (without  multi-user  interference)  and  Ph  is  the  probabil¬ 
ity  of  a  symbol  being  hit  (interferred  with  by  another  user’s  transmission,  given  by  [8] 


Using  (A. 3),  (A.2)  can  be  written  as 
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m  =0  l.  1  )  *- 

(A. 4)  has  the  form  discussed  above.  We  now  use  (A. 4)  and  the  aforementioned  exponen¬ 
tial  approximation  to  -r— —  in  (A.l)  to  yield 

i  +1 
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This  exactly  the  expression  given  as  eq.  (6)  in  Section  II. 

The  expression  for  RS(n,k)  Reed-Solomon  coding  with  erasure/error- 
correction  decoding  can  be  obtained  in  a  similar  fashion.  In  this  case,  Pc(j-i )  is 
given  by  (see,  for  example,  [9]) 


Pe(j- 0=1-  E  [J  [  „  ]  [p»  (7-*  )  ]?  [c*  (J -*  )  [!-P#  O' -*>-«•  O' -*)]  *~?’(A.6) 

e  +i<2r+>? 

where  e  —n  -k  ,  e,  ( j -t )  is  the  symbol  erasure  probability  given  by 


«,(/-*■  )=i-u-n)y'‘‘. 


and  p^O'-i),  the  symbol  error  probability,  is  given  by 


4 


for  codes  over  GF (M1)  [the  P0  would  become  1-(1-F0)m‘  for  codes  over  GF{Mm‘)\. 
Substituting  into  (A.l)  and  proceeding  as  above,  we  obtain 
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(A. 7) 

For  binary  convolutional  codes  with  Viterbi  decoding  and  hard  decisions, 


Pc(j-i )  can  be  bounded/approximated  by  (see,  for  example,  [10]) 


00 
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(A.8) 


where  rc  is  the  code  rate,  djree  is  the  free  distance  of  the  code,  Wp  is  the  total  informa¬ 
tion  weight  of  all  sequences  which  produce  paths  of  weight  p,  P0  is  as  described  above, 


p{j-i)=l-{l-P0){l-Phy  \  P=“.  and  e,  (j  -i  )=1-(1-PA  )’  f.  P2{n\q)  is  defined  as 
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It  should  be  noted  that  the  summation  over  n  in  (A.8)  generally  involves  only  a  few 
non-zero  terms. 

Substituting  (A.8)  into  (A.l)  and  proceeding  as  above,  the  following  result  is 
obtained: 
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where  for  the  case  in  which  no  side  information  is  available, 
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if  /i  is  odd,  and 
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if  /i  is  even;  whereas  for  the  case  in  which  side  information  is  available, 
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(A. 11) 

Note  that  P2(fi;—)  and  P2(/i-X;F0)  have  not  been  expanded,  since  each  is  a  constant  for 
2 

a  given  value  of  fi  or  (//-X),  respectively. 
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Appendix  B  :  Derivation  of  Expressions  for  Coded  MFR  Systems 

As  in  Appendix  A,  above,  we  will  show,  in  detail,  how  closed-form  results  have  been 
obtained  for  the  same  codes  treated  in  Appendix  A,  but  in  an  MFR  system.  The  expres¬ 
sions  for  Pc(j~i )  are  exactly  as  stated  in  Appendix  A  and  will  not  be  repeated  here. 
The  result  for  binary  convolutional  coding  with  side  information  available  at  the  decoder 
will  be  derived  in  detail,  and  the  other  results  will  be  stated,  for  the  sake  of  complete¬ 
ness. 

In  Section  III  the  following  formula  is  derived: 
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As  in  Appendix  A,  we  will  use  the  three  term  exponential  approximation  to 

3 

c  u6„.  For  convenience,  let  use  define  I„(r0,90 )  by 

i*=i 
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It  is  clear  that  this  includes  all  that  varies  from  one  code  to  another,  and  therefore,  all  of 
the  results  with  be  stated  in  terms  of  /,/ro,0o).  Note  that,  by  definition, 
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P(X^Y)  =  p[l-e  2J(l-p)ScJ0  j  Ju{roeo)fr^rQ,eo)deodrQ, 
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and,  since  expected  forward  process  is  computed  by  including  rocos0o  in  the  integrand, 

(  3  R  l 

^  P^l_e  I  (l-P  )  E  c  i'Jq  J  r  0cos90  I  v{r  0,d0)  f  r  j(r  0,60)d  d0dr  0. 
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2 

Let  us  begin  with  the  case  of  binary  convolutional  coding  with  side  informa¬ 
tion  available  at  the  decoder.  The  expression  for  Pc(j-i )  is  given  in  eq.  (A. 8).  Before 

beginning  the  actual  analysis,  let  us  first  observe  that  since  p  —  —  and  P0  is  a  constant, 

the  factors  involving  P2(;)  are  independent  of  ( j-i );  thus,  it  will  not  be  necessary  to 
expand  them  (as  opposed  to  the  case  in  which  no  side  information  is  available  to  the 
decoder;  the  expression  is  of  the  form  P2(",p(j-i))  and  must  be  expanded). 

To  start,  let  us  manipulate  Pc  ( j -i)  into  a  desirable  form. 


Pc  (i-0 


1  -  rc  E 


[eeO-,)]V2^;||+^  "j  [e»O-o]M”V2^-X;P0j 


=  1  -  r. 


li- 

+  e 

X=c 


|J“j  )'_i  ]x  [o-nr  ]“'Vs(/.-x;p0)  J 


=  1  -  rc  E  w / 

dfrce 


3 


+li(')(»)<'i,^('‘-x:n)[<i-nrX"]"'l 


(B.3) 

We  now  have  Pc(j-i )  depending  on  ( j-i )  exponentially,  as  desired. 

We  can  now  substitute  (B.3)  into  (B.2)  and  derive  I 0,60)\  however,  the  expression 
would  be  rather  long.  Instead,  we  will  treat  the  three  terms  of  (B.3)  individually  (we 
consider  the  leading  1  as  a  term).  Looking  at  the  first  term,  we  have 
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From  the  second  term,  we  get 
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Since  the  analysis  of  the  third  term  is  nearly  exactly  the  same  as  that  of  the  second 
term,  the  result  will  simply  be  stated: 
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/„(ro,0o)  is  obtained  by  taking  Ij(r0,60)  =  /„(1)  ( ro,0o )  -  /j2)  (ro,0o)  -  /j35  ( r0,60 ). 

For  the  other  three  coding  schemes  considered  here,  the  results  are  as  follows: 
1.  For  RS(n  ,k  )  Reed-Solomon  coding  with  error-correction  decoding, 
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2.  For  RS(n,k )  Reed-Solomon  coding  with  erasure/error-correction  decoding, 

use  /„(1)  (r O,0O)  in  (B.4)  and  I „(r  0,60)  =  /„(1)  (ro,0o)  -  /j2)  (ro,0o),  where 


rja)(r0.«o)=  S  E  E  I  (-i)""" 

/i=0  <  W  =o  v  ^  J  v  v  M  J 

t)=max<  0,e  +1-2 /i  > 


exp 


N 


-pN'  (ro»^o)|  1-e  2  )-[l-.P1(JR  )](1-PA 


(B.8) 

Note  that  for  both  RS  ( n  ,k )  coding  cases,  the  results  are  shown  for  codes  over  GF{M1). 
For  codes  over  GF(Mm‘),  replace  (1-P0)  in  both  cases  by  (l-P0)mc,  and,  in  (B.8), 
replace  the  factor  of  P £  by  a  factor  of  [1-(1-F0)m‘ ]p. 

3.  For  binary  convolutional  coding  with  no  side  information  available  to  the 
decoder,  we  have  I^(r0,60)  as  in  (B.4), 
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and,  if  p  is  eveni, 
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and,  as  before,  Iv{r0,60)  =  /  (r o,0o)  -  I ^  (r0,60). 
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